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o 

CN ' The motion of molecular motor is essential to the biophysical functioning 

^ . of living cells. In principle, this motion can be regraded as a multiple chemi- 

cal states process. In which, the molecular motor can jump between different 
chemical states, and in each chemical state, the motor moves forward or back- 



"^ ! ward in a corresponding potential. So, mathematically, the motion of molecular 
P-i, 

Q ' motor can be described by several coupled one- dimensional hopping models or 

i-^ . by several coupled Fokker-Planck equations. To know the basic properties of 

O I molecular motor, in this paper, we will give detailed analysis about the sim- 

*c/2 ■ 

>-,' plest cases: in which there are only two chemical states. Actually, many of 

^ ; 

Ph. the existing models, such as the flashing ratchet model, can be regarded as a 



two-state model. From the explicit expression of the mean velocity, we find 
that the mean velocity of molecular motor might be nonzero even if the po- 



(N 
> 

in 

f.^^ [ tential in each state is periodic, which means that there is no energy input to 



the molecular motor in each of the two states. At the same time, the mean 



^ . velocity might be zero even if there is energy input to the molecular motor. 

Generally, the velocity of molecular motor depends not only on the potentials 

v~j ' (or corresponding forward and backward transition rates) in the two states, 

H ! but also on the transition rates between the two chemical states, 

c^ , 
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I. INTRODUCTION 



Molecular motors are biogenic force generators acting in the nanometer range, and 
converting chemical energy into mechanical work [l|, |2] , which play essential roles in 
eukaryotic cells |3kZ|. In the super family of molecular motors 8|, the most extensively 



studied ones are conventional kinesin 



24- 



29|, and FoFi-ATPase |30 



inesin [9 



16| . cytoplasmic dynein |17H23|. myosin V 



The conventional kinesin can walk hand-over- 
hand along microtubule about 1 um to the plus end direction of the microtubule before 



its dissociation from the track 

6-8 pN [id, iii y , y, I43 



37 



with step size 8.2 nm 40l-l42| and stall force 



47j , which is independent of ATP concentration jlO| . In 



10 



saturating ATP solution, its zero load velocity is about 700—1000 nm/s 
Cytoplasmic dynein also can walk hand-over-hand along microtubule with average 
step size 8.2 nm 



22 



49 



52|, but to the minus end direction 18 1. Recent experimental 



data indicate that its stall force is also about 6—8 pN 49, 



52 



53| . and independent 



of ATP concentration 49|. To the dynein which is purified from mammalian animals 



18 



54 



55| . Myosin V is also 



its maximal velocity is also about 700—1000 nm/s 

a processive motor but walks along actin filaments with average step size 36 nm, 

and ATP independent stall force 2—3 pN 



56 



6l|. ATPase consists of two portions 



Fq and Fi connected to a 7 shaft. It can use the proton- motive force across the 
mitochondrial membranes to make ATP from ADP and Pi, and also can use ATP to 



drive the rotation of the 7 shaft 62|. Recent experiments found that there are also 
many other molecular motors that can move processively, such as kinesin CENP-E 
63|, myosin VI 64J467l|. myosin Vila 68|, myosin IXb [69|, myosin XI 70|, and T7 
DNA helicase [Til- 

There are many mathematical models to describe the motion of molecular motor, 
.uch a. Pokker-Planck equation g BQ. Langevm equation fl, and ,na=te. equa- 
tion [9|, 176H79J . However, so far, almost all of the explicit formulations of biophysical 
properties of molecular motor, such as mean velocity [2|, |80|, effective diffusion con- 
stant |8ll,l82|, and mean first passage time 83|,l8^, are obtained by employing one-sate 
models, in w hich the molecular motor moves along its track in one tilted periodic po- 
tential 108| . One of the basic properties of such models is that the mean velocity of 



molecular motor does not vanish as long as the input energy is positive. These models 
and their corresponding results are valuable to describe the tightly mechanochemi- 
cal coupled cases of motor motion. However, recent experimental data indicate the 



motion of molecular motors, including conventional kinesin |43l . ISSHSSl]. cytoplasmic 
dynein [89|, myosin II [90|, ISjJ, and Fi-ATPase |92] are usually loosely coupled to ATP 
hydrolysis, i.e., the input energy might be nonzero even if the mean velocity vanishes. 
To study these loosely coupled cases, it is necessary to use multi-state models. In 
fact, the multi-state models have been used by some authors [73|, l93|, l9J]. However, 



it is hard to get meaningful explicit results for the general A^-state models. Usually 
the numerical calculations are employed 951-197 1. 



In this paper, we will give a detailed theoretical analysis to the two-state models. 



Actually, the two-state models have most of the essential properties o: 
multi-state models, and they have been used in many studies 



95 



98 



the general 



104J . There 



are two different forms of two-state models: (1) two coupled one- dimensional hop- 
ping models, and (2) two coupled one-dimensional Fokker-Planck equations, which is 
equivalent to two coupled Langevin equations (in fact, it also can be verified that, 
any one-dimensional ho pping model can be well approximated by a one-dimensional 
Fokker-Planck equation 105| ). In the following, we will give the explicit formulation 
of the mean velocity of molecular motor by using two coupled one-dimensional hop- 
ping models and two coupled one-dimensional Fokker-Plank equations respectively. 
From this formulation, the stall force, i.e., the external load under which the mean 
velocity vanishes, can be obtained. We find that, the mean velocity, and consequently 
the stall force depend not only on potentials in the two states (or corresponding for- 
ward and backward transition rates), but also on transition rates between the two 
chemical states. In general, part of the input energy will dissipate into the environ- 
ment, and so the energy efficiency, i.e., the ratio of mechanical work done by the 
molecular motor to the input energy, might be far less than 1. For example, the 
mean velocity might be zero even if the input energy in each state is nonzero. 

The organization of this paper is as follows. In the next section, the two cou- 
pled one-dimensional hopping models are discussed, and then in Section III, the two 
coupled Fokker-Planck equations are analyzed. In each models, three special cases 



are further analyzed: (1) The motor can jump between the two chemical states at 
only one position. In fact, the properties of this special case are much similar as the 
usual one state model. At steady state, there is no energy input to the molecular 
motor during its transition between the two chemical states. (2) The motor can jump 
between the two states at two positions. This special case has the typical properties 
of the general cases. (3) One of the two potentials is constant, or all the correspond- 
ing transition rates (forward and backward) are equal to each other in one of the 
two states. This special case corresponds to the flashing ratchet model of molecular 
motors. Finally, the results are briefly summarized in Section IV. 

II. TWO COUPLED ONE-DIMENSIONAL HOPPING MODELS 

The two coupled one-dimensional hopping models are schematically depicted in 
Fig. [H In which, the forward and backward transition rates in state 1 are denoted 
hj Fn {n -^ n + 1) and Bn (n — )■ n — 1), the forward and backward transition rates in 
state 2 are denoted by /„ (n — )■ n + 1) and bn (n — )■ n — 1), and the transition rates 
between the two states at position n are denoted by w" (state 1 — )■ state 2) and u^ 
(state 2 — )■ state 1). Under the assumption of periodicity, we have 



FlN+n 


= Fn, 


BlN+n — Bn, 


JN+n 




flN+n 


^ Jni 


blN+n = bn, 


u;f+" 





where / is an integer number /, A^ is the period of hopping models. Let -P„(t) be the 
probability of finding molecular motor at position n of state 1 (denoted by 1„) at 
time t, and pn{t) be the probability of finding molecular motor at position n of state 
2 (denoted by 2„) at time t. Then the evolution of probabilities -P„(t) and p„(t) are 
governed by the following master equations: 



j/n{t) = Fn-lPn^lit) - (i^„ + 5„)P„(t) + Bn+lPn+l{t) 

-ulPn{t)+Uji~Pn{t), 
-Tj-Pnii) = fn-lPn-lit) - (/„ + &„)p„(t) + 6„+lPn+l(t) 

+ LO^nPnit) - tU%{t), U = 0, ±1, ±2, ■ ■ 



(2) 



Let 



Pnit) = J^ PiN+n{t), Pn{t) = J^ plN+n{t), 



l=—oo 



then, at steady state, -P„ and p„ satisfy [106 1 

Fn-lPn-l — {Fn + Bn)Pn + Bn+lPn+1 — ^nPn + l^nP" = 0' 
fn~lPn-l - (/n + K)pn + &n+lPn+l + ^IPn " W^^Pn = 0, 

with n = 1, 2, ■ ■ ■ , A^, and the total flux of probability. 



(3) 



(4) 



Jn+\ — {FnPn — -Bn+l-Pji+l) + {fnPn — K+lPn+l), 



(5) 



is constant, i.e. J^+i = J for n = 1,2,- ■ ■ ,N. 
From the first equation of (j4]), one sees that 



Pn — [{Fn + Bn + Uj'^)Pn — -^n-l-Pn-1 " B n+lPn+l] / ^"d ■ 



(6) 



Substituting ([6]) into ([5]), one can easily verify 



J — ^n-lPn-l + CnPn + Dn+lPn+1 + -E'n+2-Pn+25 



(7) 



where 



^n — —fn+lFn/^d ) 

a = [fn{Fn + 5„ + a;,")]/u;," + F„ + bn-,iFju^,+\ 

Dn = -[fin + /„-lfin/u;r' + ^nifn + ^n + U;:)^: 
-E'n = bn-lBn/(^d~ ■ 



By (jTj) and routine analysis, we obtain 



Pi = X,J + Y,P, + ZiP2 + W,P 



3, 



(9) 



where 



X 



1 



TV 



An 
Xn~i ~ 
Zn-i 



Yi 



N 



'An' 



D, 



^N 



An' 



W, 



E. 



N 



An' 



An-1 

Cjv-D2 

An-iA 



1-S^)> Yn^^ 



CjyCi 



£>i 



E2 



N-l^N J^N-1 



Wn-1 



An-iAm -4jv_i ' 
_ CmE'j 



(10) 



fill 



Yn-2-- 
Zn~2 - 



1 



Cn-1 



Cn-iC^ 



N 



D 



N 



An-2 An-2An-1 An^2An~iAn AiV-2^Af ' 

El Cn-iDi DnCi Cn-iCnCi 



An-2 ^Af-2^Af-l An^2An An^2An-iAn' 

Cn~iE2 , DnD2 Cn~iCnD2 



w, 



DnE^ Cn-iCnE^ 



N-2 



An~2An Aj\f_2A]\f_iAi^' 
andioT l<k< N - 3, 



Yk-- 
Zk - 



1 — (C'fe+i-^fc+i + Dk+2Xk+2 + Ek+^Xk+'i) 

Ak 
Ch+iYh+i + Du+2Yk+2 + Eu+^Yi 



^fc+S-ffc+S 



Ak 

Ck+lZk+l + Dk+2Zk+2 + Ek+3Zk+3 

A'k ' 

Ck+iWk+i + Dk+2Wk+2 + Ek+sWk+s 
Ak 



Then, for z = 1, 2, 3 in equation (j9]), we obtain the following equations 



:i2) 



(13) 



AP = JX, 



(14) 



where 



A 



( 



\ 



1 - Fi -Zi -Wi 
-Y2 I-Z2 -W2 
-Fa -Z3 l-W; 



\ 



J 



P 



/p.^ 



V^/ 



X 



X2 



(15) 



So, Pi = A^ for i = 1,2,3, with P = (A, A, A)^ satisfy AP 
Pi, ioT 3 < i < N, can be obtained by Eq. ([9]), 

Pi = (X, + yA + ZA + W-A)J =: PiJ, 

and therefore, pi can be obtained by Eq. ([6]), 



X. Consequently, 



(16) 



A =[(P, + Pi + a;:)Pi - P.-iP.-i - P,+iP,;+i]M, 
= [(P, + Bi + ul)Pi - P.-iA-i - P.+iA+i] J/c^, 



(17) 



The probability flux J in Eqs. ( !T6|) (fT7j) is determined by the normahzation condition 

N 

N 



i=l 



"""^^^a-4.V..fi-i)A. 



w^ 



^:, oj\ 



^:, oj\ 



Specially, if w^ = Ua and w^ = Ud are constants, then 

J _ i^d UJd 



n. 



N N 



(18) 



(19) 



A. Special case I: w^ = Wj = for 1 < i < A^ — 1 



For convenience, we denote uj^,uj^ by uja,ujd respectively (see Fig. |2]). For this 
special case, the steady state probabilities P„, p„ satisfy 



FnPn — BiPi — FiPi — B2P2 — ■ ■ ■ — Fn-iPn-i — BnPn —'■ J, 
JnPn - hPi = fiPi - &2P2 = ■ ■ ■ = Jn-iPn-i - bNPN =■ j- 
It can be readily verified that 



(20) 



Pk — XkPN — YkJ, 



(21) 



with 



Specially, 



which implies 



k k 






N 



F^-i 



1 A^K. 



N 



n^K-^En^h'. 



,i=l 



1=1 j=i ■' 



N N 



F„ Z^ 11 R, 



Fn 



N 



i=lj=i 



N 

n 



-J. 



n%^-i 



(22) 



(23) 



(24) 



Combining (EI]) (EJ]) and ([21]), one finds 

/ / N N 

^^ Enf i(n 



Pk 



Fn 



i=lj=i 



i=l 



F\-j_ 

Br 



\ 



\ 



N 

n 



N N ^ 

Enf 

8=1 j=l 



F\-i_ 

Br 



n 



-1 



1 k k J-, 

1=1 j=i ■' 



J, 



k k 

En 

1=1 j=i 



B, 



N 



N 

n 



n%r-i 



(25) 



Using the periodic conditions ([1]), one can verify that 

N+k N+k „ 



t E n 



Ffc 



i=k+l j=i 



N 



-J. 



nft-1 



j=i 



Using the same method, the probabihty pk can be obtained 



Pk 



N+k N+k 

i E n 

i=k+l j=i 

iv 



(26) 



(27) 



At steady state, Wa-P/v = <^dPN, which imphes 



Therefore, 



N N 



Fn 



i=lj=i 



B, 



N 

n 

i=l 



B, 



COn 



N N 

srEn 

1=1 J=l 



N 



nt-i' "" 



i=l 



J =: EJ. 



N+k N+k 

i E n^ 



N N 



fk . 



j£q_ Y^ TT £i 

F« ^ 11 B, 



Pk 



i=k+l j=i 



Fn 



N 

n 



-GJ, where O 



j=l jr=i 



nf-1 



N N ^ 

^Enlf 



fN ^ 



i=lj=i 



N 



(28) 



(29) 



Since Pk, Pk satisfy E (-^fe + Pk) = 1, from (l26l [291) . the probabihty fiux J can be 

fc+i 
obtained as follows 

N 



J 



ni-1 



i=l 



(30) 



0^ + ^V' 



where 



N N 



^d Y^ fr fj 
fN ^ -li 6, 



$ 



N N 

1=1 J=J 



El 



Af / 7V+A: N+k . 

fc=i y i=k+i j=i ^ 

N ( N+k N+k 



(31) 



fc=i 



i=fc+l j=j 



So the total flux of this system is 
J + j = (l + S)J = 



TV \ 



1 + ^ e 



V 



N 

1=1 



J 



N 

n 

i=l 



B, 



N 

^j=i 



0\I^ + $v^ 



/ 



(32) 



Combining ( l26l) ( l29l) and ( l30l) . the probabilities Pfc and p^ can be obtained as follows 

N+k N+k „ -, _ Af+fc N+k 



0^ + $^ Ffc 



^En 



i=k+l j=i 



N 






Pk 



$ 



0^ + ^v- A 



^En 



i=fc+l j=j 






(33) 



By fl32|) . one easily sees that, if H {^ 



N 



fl ^ = 1, then J = j = 0, consequently 



the total probability flux J + j = 0. In other words, for this special case, if there is 
no energy input to the molecular motor in each state, then the mean velocity would 
be zero. But the reverse does not hold. Note, the potential changes in one period 
of state 1 and stat e 2 are AGi = A;sTln (ni=i ^) ^^d AG'2 = kBT\i\ (nj=i b' 
respectively J9l. Il07 |. 



10 



B. Special case II: ui = uj'^ = for i j^ M,N 



For convenience, we denote u^ , u^ by fla, ^d, and u^ , ujf by Ua, cOd (see Fig. |3]). 
At steady state, Pk,Pk satisfy 

FnPn — BiPi = FiPi — -B2-P2 = ■ ■ ■ = Fm-iPm-i — B^Pm ='■ Ji, 
FmPm — Bm+iPm+i = ■ ■ ■ = Fn-iPn-1 — B^P^ =: J2, 
InPN - biPl = flPl - &2P2 = ■ ■ ■ = JM^lPM-l - buPM =■ ji, 
ImPm — bM+lPM+l = ■ ■ ■ = Jn-iPn-i — ^nPn ='■ J2, 

UaPu + ^aPN = ^dpM + ^dpN, ^ ' 

J2 = Jl + ^a-P/V — ^dPN, 

J2 = jl - ^aPN + ^dPN, 

N 

Y,{Pk+Pk) = i. 

k=l 

From the first equation in fl5^ . one can easily get 
Pk- 



^ T7< \ 1 / ^ ^ 17 



\i=l 



Bi 



At the same time, from the second equation in (j3 

k 



P,= 



F\_i 
B, 



n ^-^\pm 



Vi=M+l 
\i=M+l '■ 



i=M+l j=i •' 



i^.-l 






vi=l 



B^ 



(k k \ 

i=A/+l j=i ^ J 



\i=Af+l j=t •' 



.4 = 1 

'Fk 



Pn 

k k 



= :{Rk- naSk)PN + ^dSkPN - TkJl, ioiM+l<k<N. 
In particular, P/v = {Rn — ^aSN)PN + ^dSNPN — T^Ji, which gives 

{Rn — ^aSN — 1)-P/V + ^dSNPN 
J^ — 



T, 



(35) 



(36) 



N 



(37) 



Substituting (1371) into (l35l) (l36l) . we obtain 



11 



-Pfc — GuPn + HkpN, 



(3J 



where 



Gi 



-f?^fc 



i?fc - §(i?7V - ^aSN - 1), 1 < A; < M, 

i?fc - naSk -^{Rn- fia^TV - 1), M + 1 < A; < iV, 

-^fid^jv, 1 < A; < M, 



M+l<k<N. 



(39) 



Similarly, 



Pfc = gkpN + /ifc-Pj 



Af; 



(40) 



where g^, h^, and the corresponding r^, s^, t^ in expressions of Qk, h^, can be obtained 
by replacing Fj, Bj, Qa, ^d in the expressions of Rk, Sk, Tk, Gk, Hk with fj, bj, Qd, ^a 
respectively. Combining fl55]) fffiJ]) and the fifth equality in flM|) . we have 



UJa{GMPN + HmPn) + ^a-PAf = ^d{gMPN + /iM-P^v) + f^dPAf, 



(41) 



I.e. 



(fia + CUaCAf " UJdhM)PN = (^d + ^dS'A/ " UJaHM)p 



N- 



(42) 



So 



^a + ^aGu — ^dhjvl J-, U 

ild + UJdQM - ^aiiM V 



From ([MD (HOD and (gSD, one finds 



(43) 



Pk — \Gk+ Y^k 1 Pn, Pk — {hk + —gk ] Pn- 



(44) 



In view of the last equation in (l3ll) . one gets 



J2U^k + hk) + -{gk + Hk)\ 

k=l ^ ^ 



PN = h 



(45) 



which implies 



P 



N 



N 

Y. 

k=l 



E [{Gk + hk) + ^{gk + Hk 



(46) 
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BydSZD (E 


m^ 


m. 


we have 














Ji = 


(-Rat - 


■ ^aSN — 


1) + ridSNy 




N 

TnEH 
k=l 


Gk + h) 


+ M,(^g, + Hk)] 


Similarly, 


one 


can 


verify that 












Ji = 


{rN- 


^dSN - 


l)^ + QaSN 




N 

tN E [(< 
k=l 


S'fc + h) 


+ ^{g, + H,)] 



(47) 



(4J 



Therefore, the total flux of this special case is 

Ji+3i — J^ • (49 j 

TivtjvE [{Gk + hk) + '^{gk + Hk)\ 

k=l 

More specially, if Rn = tat = 1, then the total probability flux is 

fc=l 



^a^dfM — ^a^dRu (Sn S^ 



m_ 



E[(G'fe + /ife)V^+(^fc + ^fc)f/] 



(50) 



fc=i 



where 



(51) 



t/ =Qa + ^aRM + ^a^aSNTM/TN + Ud^aSNtM /^N > 0, 
V^ =lld + WrfrM + UdfldSNtM/tN + OJafldSNTM/TN > 0. 

So the direction of probability flux is determined by the sign of ( f^ ~ f^ ) and 
{Qai^drM — ^S^dRhi)- One can see that, Rn = r^r = 1, i.e., AGi = /S.G2 = 0, does 
not read the mean velocity vanishes. 

To better understand the inter-state transition rates dependence of the total prob- 
ability flux, we assume that 

{^a, ^d, <^a, ^d) = ^{^a, ^d, ^a, ^d)- 

It can be verified that the total probability flux J := Ji + ji in fl50|) increases mono- 
tonically with parameter A. If A = then J = 0. If A —t- 00, them J tends to 

^a^d^M — UJoVLciRm f Sn_ _ ^N_\ /p-2^ 
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where 



N 

fc=l 



-)Rk + 7^ \ : )rk 



t 



N 



Tn 

N 



E 



+ {^a'^d'^M — ^a^dRpI 



.k=l 

N 



TkSN 



+ {VLaUdTM - UJa^dRu) 

I : \ 7^ 1 2_^{VtdRk + ytaTk 

N 



^ ) + E (^fc - Sk 



tN J 



fc=M+l 



k=\ 



TkSN _ tkSN\ V^ 



tN J 



[sk — Sk 



k=M+l 



C. Special case III: w^ = w^ = for i / M, N, and fi = bi = f ior I < i < N 



As pointed out in the Introduction, the flashing ratchet model can be regarded as 
one example of this special case. For this more special case, we have 



k 
Tk = 1, tk = -, for 1 < /c < A^, 

and Sk = {k — M)/ f for M + 1 < A; < A^. It can be easily verified that 



9k 



. (N-M)k ^ 
l, M(N-k) ^ 

(Af-M)fc n 
M{N-k) Q 



A^/ 



1 < fc < M, 

M + 1 < A; < AT, 

1 <k< M, 

M + l<k<N. 



So 



f/ =(]„ + CJaGM - ^dh 



M 



V -LN J-N J ^V/ 

V =rid + UJdgM — ^aHM 

=ild + CJd I IH TT-: ild I + —SNUjJld- 

Moreover, if Rn = 1 then the total probability flux is 

Qa^d - ^a^dRM f N - M _ Sn_ 

N Tn 



Ji+Ji 



(53) 



(54) 



(55) 



(56) 
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where 



M 



M 



N 



A = ^ {VRk + Urk) + {VtaUdTM - UJa^dR 
fc=l 

M 

+ MU + V^Rk. 



M) 



EfTkSj^f tkS]\f\ Y^ / o 

k=l ^ ^^ ^^ ^ fc=M+l 



^^ fc=l k=M+l 



2Nf 



k=l 



III. TWO COUPLED FOKKER-PLANCK EQUATIONS 



The general two coupled Fokker-Planck equations are as follows 



dtP =Dd^i/3VlP + d^P) + ujd{x)P - Ua{x)P, 
dtp =Dd^{(5V2p + 9^p) - ujd{,x)p + uja{x)P, 



oo < a; < +CXD, 



(57) 



where D is free diffusion constant, (5 = l/ksT with ks is Boltzmann constant, and 

T is absolute temperature, P{x, t) and p(x, t) are probability densities of finding 

molecular motor at position x at time t and in states 1 and 2 respectively, Vi, V2 are 

(tilted) periodic potentials with period L. uja{x\ujd{,x) are transition rates between 

states 1 and 2 at position x |l09| . Similar as in 82|, let 

+00 +00 

P{x,t)= ^ P{x + kL,t), p{x,t)= ^ p{x + kL,t), (58) 

k=—oo k=—oo 

then P{x,t), p{x,t) satisfy 



< X < L. 



(59) 



9iP =Dd,{f3V;P + d,P) + Ud{x)p - u;,(x)P, 
(9ip =Dd^{(3V^p + 9^.p) - a;d(a;)p + oja{x)P, 

The steady state solution of ( l59l) can be obtained under the following constraints: 
P{0) = P{L), p(0)=p(L), / {P + p)dx = l, / Udpdx= UaPdx. (60) 
The corresponding probability flux is 



J=-D {i3v;p + d,P + /^^V + d.p) 



(61) 
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and the mean velocity of molecular motor is V^ = /^ Jdx = —[3D J^ (V7P + 1^'p)(ix = 
JL. If cOa{x),cod{x) are constants, Eq. ( l59l) had been discussed by Y.-D. Chen [95|] 



and it can 
method 



DC solved numerically using the similar method as the one used in WPE 



96 



97|. 



A. Special case I: uja{x) = ujd{x) = for < x < L 

For this special case, the steady state probability densities P{x),p{x) of finding 
molecular motor at position x are governed by the following equations 

Dd,{pv;p + d,p) = o, 

0<x<L. (62) 

Dd,{l3V^p + d^p) = 0, 

Meanwhile, P{x),p{x) satisfy the following boundary conditions and normalization 
constraint: 

P(0) = P(L), p(0)=p(L), LOaPiO)=tUdp{0), f iP + p)dx = l, (63) 

Jo 

where Ua = uJa{L),ujd = <^d{L). The probability fluxes in the two states are 

J = -D{PV;P + d,P), J = -D{/3V^p + d,p). (64) 

So Eqs. ( J62l) can be reformulated as 

PVlP + d,P = -J/D, 

0<x <L. (65) 

(3V^p + d,p = -j/D, 
The general solutions of (!65|) are 

P(a;) = fl re^y^(y)dy + C^ e-^^^^^\ p{x) = (-^ r e^^^^y^dy + C^ e''^^^^^), 

where the constants Ci,C2 can be determined by the periodic boundary conditions 
P(0) = P(L),p(0) = p(L): 

_ j,J^\^y^iy)dy _ lj^^e^v^(y^dy 

with AVi = Vi{0) - Vi{L). Therefore 

J f^+L pl3[Viiy)-Vi{x)]j j_ f^+I^ pl3[V2iy)-V2{x)]^ 
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From LOaP{0) = Wrfp(O), one sees 



^a ; _«Ai/, — ^d- 



1 _ e-^^^i 



1 _ e-/3AV2 



so 






(67) 



From fl66|) fl67|) and the normalization condition J^ (P + p)(ia; = 1, one can easily get 



J 






{6t 



where 



-k =Ude 



/3Vi(0) 



+ WaC 



/3V2(0) 



,/3V2(.)^, 



./3Vife) 



-l3Vi{x) 



uo 



x+L 



'i^^^^yUy]dx 



dy 



-fSViix) 



IJO 



X+L 



z^^^^y^dy 1 dx 



It can be easily found that, for this special case, the total probability flux J + j = 
if potentials V^i, V2 are periodic, i.e., AVi = AV2 = 0. From fl32|) and fl68l) . one sees 
that, the properties of this special case are simil ar as those of the special case I of 



the two coupled one-dimensional hopping models 



105J 



B. Special case II: 0Ja{x) = 0Jd{x) = for x ^ a, L 



For this special case, the governing equations of the steady state probability den- 
sities P{x),p{x) are 



Dd,{pV^p, + d,pi) 
Dd,{f3VlP2 + d,P2) 
Dd,{f3Vip2 + d,P2) 



0, 



< a; < a, 



a < X < L, 



(69) 
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with the following constraints 

PM = P2{L), Pi{a) = P2{a), 

piiO)=p2iL), pi(a)=p2(a), 

Ji = J2 + uJaPia) - Udpia), 

ji = J2 - ^aP{a) + u;dp{a), 

UaP{a) + QaP{L) = ujdp{a) + ^dp{L), 

pa pL 

/ {Pi+pi)dx+l {P2 + p2)dx = l, 

Jo Ja 

where Ua = Ua{a),Ud = uJd{a), Qa = (^a{L), Qd = ^d{,L), and 

J, = -D{I3V[P, + d^Pi) ji = D- {l3Vipi + d^pi), for z = 1, 2, 

are probability fluxes in the two states. 

The general solutions of (l69l) can be written as follows 



Fi(x) = ;^e-^^i(") /; e^^^'^y^dy, Gi{x) = g-'^^i^"), < x < a, 
F^{x) = ^e-^^i(^) XT e'^^'^^'^dy, G2{x) = e-^^^^''\ a<x<L 
f^{x) = ;^e-^^2(^) /o" e^^^^y^y, gi{x) = e-'^^^W^ < x < a, 



(70) 



Pi{x) = -Fi{x)Ji + Gi{x)Ci, pi{x) = -fi{x)ji + gi{x)ci, i = 1,2, (71) 
where 



(72) 



From (170!) (ITTj) . one can verify that Jj, jj and Cj, Cj satisfy the following equations 

- Fi{a)Ji + Gi{a)Ci = ^2(0)^2, 
5'i(0)ci = -f2iL)J2 + g2iL)c2, 

- /i(a)ji + 9iia)ci = 5-2 (a) C2, 

Ji = J2 + a;aG'2(a)C2 - Wd5'2(a)c2, 
^1 + ji = J2+J2, 



(73) 



WaG'2(a)C2 + fi„Gi(0)Ci = iOdg{a)c2 + lld5'i(0)ci, 

Fidx 1 Ji + ( / Girfx 1 Ci - ( / F2dx ) J2 



6*2 c^a; I C2 



fidx ji 



gidx ci - 



f2dx J2 



g^dx C2 = 1. 



For the sake of convenience, we rewrite equations in (!73|) as AX = B, with X 
( Ji, Ci, J2, C2, Ji, ci, j2, C2)^, 5 = (0, 0, 0, 0, 0, 0, 0, 1)^ and 



A 






Gi(0) 


F2{L) 


-G2{L) 











-F,{a) 


Gi(a) 





-G2{a) 


























^i(O) f2{L) 


-92{L) 














-/l(«) 


(71 (a) 


-92ia) 


1 





-1 


-UJaG2{a) 








^d92{a) 


1 





-1 





1 


-1 





-/Fi 


IG, 


-IF2 


IG2 


-//l 


I91 -If2 


I92 


ra 


VT 1 T 1 


fT 


pL TT 1 r 


TT 


T-i j^ r A 


J 1 1 



in which IHi = j^ Hidx, IH2 = J^ H2dx for H = F,G, /, g. Although it can be 
obtained explicitly, the solution of AX = B is very complex. So, for simplicity, we 
only discuss the special cases in which potentials Vi, V2 satisfy AV^i = AV2 = 0. By 
routine analysis, one can obtain 

Ji = F2{L)G,{a)[-uJaQdGi{a)g,{0)f2{L)g,{a)-UandGi{a)f\{a)g^{0f 
QaOo,G^iO)g,{a)y2iL) + Qai^aGmfiia)9i{0)9i{a)]/ det{A), 

J2 = -F^{a)G^i0)[-UandGlia)g^{0)f2{L)g^{a)-UaQdGlia)f^{a)g,{0y 
QaOOdGm9i{a)''f2{L) + naUJdGmf^ia)g^{0)g^{a)]/ det{A), 
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ji = -/2(L)(?i(a)[-a;„firfFi(a)Gi(0)Gi(a)(7i(0) -a;AF2(L)Gi(a)V(0) 
+naOOdFi{a)Gi{Oygi{a) + naUJdGi{0)F2{L)Gi{a)gi{a)]/ det{A), 

J2 = ^{a)g,{0)[-UandFl{a)G^{0)G,{a)g,{0)-^a^dF2{L)G^{ayg^iO) 
+QaU^,F^{a)G^{Ofg^ia) + Qai^,G^{0)F,{L)G^ia)g^{a)]/det{A), 
where det(y4) is the determinant of matrix A, and it can be proved that det(yl) < 0. 
So the total probabihty flux is 



Ji 



h =J2+J2 

= [naCUdGiiO)gi{a) - a;,fidGi(a)(7i(0)] 

X [/i(a)(7i(0)F2(L)Gi(a) -Fi(a)Gi(0)/2(L)(7i(a)]/det(A) 



[(7i(0)Gi(0)]Vi(a)Fi(a) 
det(A) 



^a^d 77TT ^ ^a^d- 



gm 



'Gi(O) 



X 



F2(L)Gi(a) /2(L)^i(a) 



(74) 



LFi(a)Gi(0) fi{a)g,{0)\ 
[<7i(0)Gi(0)]2/i(a)Fi(a 



X 



det(A) 

J a Ja 



[n,a;rfe^(^^(°)-^^('^« - a;«(],e«^i(°)-^^('^))] 



^ f.liV2(y) 



j^epy^^y) j^e^y^^y) 



Obviously, Ji + Ji > if and only if [f^^w^e'^^^^'^o^-^^Ca)) _ cu^fi^e^^^i^o^-^^^'^))] x 



J'^e/9Vi(!;) ^^ ePV2(y) 

pgfiViiy) ~ pgflV2{y) 



< 0. In view of the expression in f l50|l . one can find that, 
the properties of this special case are similar as those of the special case II of the cou- 
pled one- dimensional hopping models. The mean velocity of molecular motor might 
not be zero even if there is no energy input in each state. For this special case, the 
energy for motor m otion comes from the processes that drive the motor from one 
state to another |98l - ll01 |. 
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C. Special case III: C0a{x) = (^^{x) = for x y^ a,L, and V2{x) is constant 

For this special case, the governing equations of steady state probabihty densities 
P(x),p(x) are as follows 

Dd,{l3V;Pi + d,Pi)=0, 

< X < a, 

(75) 

Dd,iPVlP2 + d,P2)=0, 

a < X < L. 

Ddlp2 = 0, 

Its general solutions are ( 17T1) but with fi{x) = x/D,gi{x) = 1. The solution which 
satisfies the constraints ( 1701) is as follows 

Ji = -2[uJdGi{0)G2{a)DL - iOdGi{a)G2{L)DL + ndGi{0)G2{a)DL 
-ndG^{a)G2{L)DL + fida;^G'i(0)G2(a)aL - fi^WdCi (0)^2(0)02 
-ndUJdGi{a)G2{L)aL + ndUJdGi{a)G2{L)a'^ 
-UandGiia)F2iL)G2ia)DL + natUdGi{0)F2{L)G2ia)LD]/detiA), 

J2 = -2[iUandFiia)GiiO)LDG2ia) - fi,u;rfFi(a)Gi(0)LDG'2(L) 
+UdGi{0)G2{a)DL - cudGi{a)G2{L)DL + ndGi{0)G2{a)DL 
-QdGi{a)G2{L)DL + QdUJdGiiO)G2ia)aL 
-ndOOdGi{a)G2{L)aL - firfa;,Gi(0)G2(a)a' + ndiOdGi{a)G2{L)a^]/ det{A), 

ji = 2(L - a)[-ujandFi{a)G,{0)G2{a) - uJandG,{a)F2{L)G2{a) 
+naUJdF,{a)GiiO)G2iL) + naUJdGiiO)F2iL)G2ia)]D/ det{A), 

j^ = -2a[-uJandF,{a)G,{0)G2{a)-iOandGi{a)F2{L)G2{a) 

+naUdF^{a)GiiO)G2iL) + naUdGiiO)F2iL)G2ia)]D/ det{A). 

So the total probability flux is 

■Ji +jl =J2+J2 

=2{[Gi{a)G2{L) - Gi{0)G2ia)][ujdDL + QdDL + Ud^diL - a)a] 
+ 2aDF2{L)G2{a)[uja^dGi{a) - ^aUJdGiiQ)] 
+ 2D{L - a)F^{a)GM[^aUJdG2{L) - a;AG2(a)]}/ det(A). 
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More specially, if potential Vi{x) is periodic and continuous at a, then Gi{0) 
G2{L),Gi{a) = G2{a). So 

2D 

Ji+Ji =d^^[^'^^'^^i^") " ^aOOdGim[aF2{L)G2{a) - {L - a)F,{a)G,{0)] 



det(A) 

X (a r e^^^^^^dx -{L-a) T e^^^^^^dx 
2D{Gi{Q)fG2{a 



[a;afi.e'^(^^(°)-^^('^))-f^aC..l 



det(A) 



det(A) 



(77) 
Therefore, the total probability flux Ji + ji > if and only if [w^fide'^^^i^^^-^i^'')) - 
^a^d] (a/o^e^^i(^)da;-L/g''e^^i(^)da;) < 0. Similar as before, from ([56D and ([77D 



one can find the similarity between them 105 1. 



To better understand the properties of our model, we discuss the direction of prob- 
ability fluxes here. For the special case in which there are only two locations at which 
the inter-state transition rates are nonzero, i.e., the special case II, there are alto- 
gether 18 different types of probability flux. Since the states 1 and 2 are temporally 
symmetric, we restrict our discussion only on the cases in which UaP{a) —ujdp{a) > 
for the continuous model, or UaPM — <^dPM > for the hopping model. Then, there 
are altogether 9 different types of probability flux (see Fig. H]). Furthermore, if 
AVi = AV2 = 0, then there is only one type (see the figure (2, 2) in Fig. H]). On the 
other hand, if the potential V2 is constant (or fi = bi = f for the hopping model), 
then there are altogether 3 different types of probability flux (see the second column 
in Fig. H]). 
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IV. CONCLUSIONS 

In conclusion, two chemical states models of molecular motor are discussed in this 
paper. For some special cases, explicit expressions of mean velocity are obtained. We 
find that the mean velocity of molecular motor might not be zero even if both of the 
potentials are periodic, which means there is no energy input to the molecular motor 
in each of the chemical states. The energy for the motion molecular motor motion 
comes from the processes that drive the motor from one state to another. For motor 
proteins, these processes are ATP hydrolysis. At the same time, from the expression 
of mean velocity, we find that the velocity of molecular motor might be zero even if 



there exists nonzero input energy. Which imp 



usually loosely coupled to ATP hydrolysis [43 



ies that the motion of motor protein is 
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m. 
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with an effective potential. Theoretically, this effective potential can 
weighted average of potentials in each state of the multi-state model [£ 
[109] For motor proteins, u}a{x),ujd{x) depend on the standard chemical potentials and 
concentrations of ATP, ADP and ionic phosphate Pi J2|. llOO ]. 
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FIG. 1: Schematic depiction of two coupled one-dimensional hopping models. In which the 
forward and backward transition rates of molecular motor in state 1 are denoted by Fn and 
Bm and are denoted by /„ and 6„ for molecular motor in state 2, here 1 < n < A^ with 
N is the period of the hopping models. The inter-state transition rates at position n are 
denoted by oj" (states 1^-2) and uj'^ (states 2— >1). For motor proteins, w", uj'^ depend on 
the chemical potentials and concentrations of ATP and ADP. 
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FIG. 2: Special case I of two coupled one-dimensional hopping models. In which w" 



OJa 



for 1 < n < A^ — 1, and uj 



N 



Wa, ^d 



(jjd- For this special case, the mean velocity 



of molecular motor would be zero if there is no energy input in each of the two states, i.e., 
AGi = AG2 = 0. In fact, at steady state, there is also no energy input during the process 
that drives the motor from one state to another, since OJaPN = ^dPN- 
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FIG. 3: Special case II of two coupled one-dimensional hopping models. In which uj'^ = 
u}2 = ior n ^ M,N, and lo^^ = oja, OJ^^ = ujd, ^a = ^a, ^d ~ ^d- For this special case, 
the mean velocity of molecular motors might not be zero even if there is no energy input 
in each of the two states, i.e., AGi = AG2 = 0. Since there usually exists energy input 
to molecular motor during its jump from one state to another unless oJaPkl = ^dPM and 
^a-P/v = ^dPN- 
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FIG. 4: Different types of probability flux for special case II. From these figures, one can see 
that the motion of molecular motors is usually loosely coupled to the energy input process, 
i.e., there might exist energy input but without directed macroscopic mechanical motion. 
Part of the input energy will be consumed during substep oscillation. 



